home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 4 / Apprentice-Release4.iso / Languages / RLaB 1.18c / toolbox / faxis.r < prev    next >
Encoding:
Text File  |  1994-10-07  |  1.0 KB  |  46 lines  |  [TEXT/ttxt]
  1. //---------------------------------------------------------------------------
  2. //  faxis.r
  3.  
  4. //  Syntax:    faxis ( X )
  5. //        faxis ( X , T )
  6. //        faxis ( X , T, axis_type )
  7.  
  8. //  Description:
  9.  
  10. //  Faxis generates a frequency axis for FFT plots.
  11. //
  12. //  X = FFT data
  13. //  T = sampling period (optional argument)
  14. //  axis_type = type of axis to create:
  15. //              1 = Digital Rad/s [0,2pi], 2 = Analog Radians/s
  16. //              3 = Analog Hertz           4 = Normalized frequency [0,2]
  17. //
  18. //  Defaults if not specified:   T = 1, axis_type = 1
  19. //
  20. //---------------------------------------------------------------------------
  21.  
  22. faxis = function ( X, _T, axis_type )
  23. {
  24.   global (pi)
  25.  
  26.   if (!exist (axis_type)) { axis_type = 1; }
  27.   if (!exist (_T)) { T = 1; axis_type = 1; else T = _T; } 
  28.  
  29.   N = length (X);
  30.   a = (0:N-1)/N;
  31.   a = reshape (a, X.nr, X.nc);
  32.  
  33.   if (axis_type == 3)
  34.   {
  35.     a = a/T;
  36.   else if (axis_type == 2) {
  37.     a = a * (2*pi/T);
  38.   else if (axis_type == 4) {
  39.     a = 2 * a;
  40.   else
  41.     a = a * (2*pi);
  42.   }}}
  43.  
  44.   return a
  45. };
  46.